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Abstract 

We study polylithiated molecules as building blocks for hydrogen storage materials, using 
first-principles calculations. CLi4 and OLi2 bind 12 and 10 hydrogen molecules, respectively, 
with an average binding energy of 0. 10 and 0. 13 eV, leading to gravimetric densities of 37.8 and 
40.3 weight % H. Bonding between Li and C or O is strongly polar and H2 molecules attach to 
the partially charged Li atoms without dissociating, which is favorable for (de)hydrogenation 
kinetics. CLi^ and OLi^ molecules can be chemically bonded to graphene sheets to hinder 
the aggregation of such molecules. In particular B or Be doped graphene strongly bind the 
molecules without seriously affecting the hydrogen binding energy. It still leads to a hydrogen 
storage capacity in the range 5-8.5 wt.% H. 

^Computational Materials Science, Faculty of Science and Technology and MESA+ Research Institute, University 
ofTwente, P.O. Box 217, 7500 AE Enschede, The Netherlands. 

^Electronic Structure of Materials, Institute for Molecules and Materials, Faculty of Science, Radboud University 
Nijmegen, Heyendaalseweg 135, 6525 AJ Nijmegen, The Netherlands. 
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In recent years we observe a widespread effort aimed at finding alternatives to fossil fuels. The 
goal is to apply renewable resources that meet the increasing global demand for energy, while at 
the same time cut down CO2 emissions.^ Hydrogen is a potential energy carrier, since it can be 
produced sustainably and is one of the most abundant elements.^ Burning hydrogen does not af- 
fect the environment as only water is produced. On board utilization of hydrogen on a large scale 
is hampered, however, by the lack of efficient hydrogen storage systems. Some of the methods 
suggested so far are high pressure gas tanks, low temperature liquid hydrogen tanks, chemical 
solid state storage such as metal-hydrides, and physisorption of hydrogen by materials with a high 
surface area.^^ Physisorption offers the possibility of storing hydrogen in molecular form. From 
kinetic considerations this is advantageous over chemical storage in atomic form, which requires 
dissociation of the hydrogen bond and the formation of a hydride. Materials that are mainly based 
upon carbon are cheap, lightweight, and can have large internal surfaces. The bottleneck of ph- 
ysisorption by such materials lies in the weak interaction of hydrogen molecules with the host 
materials.^ As a result, very low temperatures or high pressures are necessary to attain meaningful 
storage capacities."^ One can imagine encapsulating hydrogen in carbon nanocages, but the appli- 
cability of such systems is unknown at present.*^ At the moment none of the above methods have 
achieved all requirements of a hydrogen storage system. 

Previous studies have shown that isolated metal atoms or ions can bind hydrogen molecules, 
where the type of bonding depends on the metallic species. Transition metal atoms bind hydrogen 
molecules through weak multicenter chemical bonds called Kubas bonds. I^^^ Alkali and alka- 
line earth ions bind hydrogen molecules through electrostatic (charge-multipole) and polarization 
(charge-induced multipole) interactions. Of particular interest is the lightweight Li+ ion, which 
can bind up to six hydrogen molecules with an average binding energy of 0.19 eV/H2.'^ This 
is close to optimal, as the binding energy should be around 0.15 eV/H2 for storage of hydrogen 
gas at moderate pressures.^ For this type of bonding it is essential that Li is ionized. More- 
over, Li atoms need to be prevented from clustering, as, for instance, steric hindrance would 
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result in a significant reduction of the amount of hydrogen molecules that can be bonded. Re- 
cent research therefore focuses on the hydrogen storage properties of Li-doped carbon nanostruc- 
turesJ ^^'^^'^^'^^'^^ ' ^^'^^'^^'^Q'^^ ' ^^'^^'^^'^^ 'A hmited amount of Li can be introduced by doping without 
the Li atoms clustering, consequently restricting the wt.% of adsorbed hydrogen significantly. 

Polylithiated molecules are species that contain a large density of Li atoms. Such molecules 
have been studied experimentally and theoretically for a number of years, but never with an eye on 
their potential for hydrogen storage. The main focus has been on the interaction of single atoms 
such as C or O, with multiple Li atoms. The first polylithiated molecule, OLis, has been identi- 
fied by mass spectroscopy by Wu et al!^ and its bonding characterized quantum chemically by 
Schleyer et al.^ The latter group has also predicted the existence of the carbon based polylithi- 
ated molecules CLi5 and CLi6,l^ which has later been confirmed experimentally by Kudo.^^ These 
studies have stimulated extensive research into the properties of CLi^ and OLi^ molecules. l-^^ * -^^*^^ * -^-^*-^^*^^ * ^^ ! 
C— Li and O— Li bonds are strongly polar, resulting in a significant positive charge on the Li atoms. 
This should facilitate electrostatic bonding to hydrogen molecules, as discussed above. 

In this paper we study the hydrogen storage properties of the polylithiated molecules CLi^, 
n = 3-5 and OLi^, m = 1-4. We find that each Li atom can bind multiple hydrogen molecules 
and that the binding energy increases with the partial charge on the Li atom. In particular, a CLi4 
molecule binds up to 12 and an OLi2 molecule up to 10 hydrogen molecules, with binding ener- 
gies in the range of 0.10-0.14 and 0.13-0.21 eV/H2, respectively. The corresponding gravimetric 
storage densities are 37.8 wt.% H for CLi4H24 and 40.3 wt.% H for OLi2H2o. Such extremely 
high densities make polylithiated molecules fascinating building blocks for hydrogen storage ma- 
terials. Aggregation of CLi^ or OLi^ molecules would be detrimental for the reversibility of hy- 
drogenation. Therefore we explore the possibilities of attaching polylithiated molecules to carbon 
nanostructures, in particular graphene, in order to immobilize them. Using dopants such as B or 
Be, one can increase the binding of the molecules to the substrate, without weakening the binding 
of hydrogen to the molecules. This scenario leads to several interesting structures with gravimetric 
densities > 5 wt.% H. 
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Computational Details 

The calculations are performed at the level of density functional theory (DFT), using the PW91 
generalized gradient approximation (GGA),'^ and the projector augmented wave (PAW) method,^^'^ 
as implemented in the Vienna ab initio simulation package (VASP).'^^'^The PAW potentials treat 
the H l^", Li and Be l^", 2s, B, C and O 2s, 2p as valence electrons. With a plane wave kinetic 
energy cutoff of 518 eV the total energies are converged on a scale of 10~^ eV/atom. Geometries 
are optimized without any symmetry constraints, using the conjugate gradient (CG) algorithm. 
Convergence criteria of 10~^ eV and 0.02 eV/A are applied for the energy difference between 
consecutive electronic steps and the total remaining forces on the atoms, respectively. Single 
molecules are treated within a cubic box with a cell parameter of 20 A, using periodic boundary 
conditions and F-point sampling. In the graphene systems to be discussed below we use a 4x4 in- 
plane graphene supercell and a 3 x3 regular k-point mesh to sample the Brillouin Zone (BZ). The 
final total energies and charge densities of the optimized structures are obtained using the corrected 
linear tetrahedron method.!^ 

Results and Discussion 
Polylithiated Molecules 

First we analyze the interaction of molecular hydrogen with CLi„ and OLi^ molecules. As the 
electron affinities of C and O are much larger than that of Li, the C— Li and O— Li bonds are very 
polar. If C and O become hypercoordinated by Li atoms, i.e. CLi^ with > 5, or OLi^ with m > 3, 
the Li coordination shells become increasingly metallic and the average charge on the Li atoms 
decreases. The binding of hydrogen to these atoms is dominated by electrostatic and polarization 
interactions. Therefore, to maximize hydrogen binding energies, the charge on the Li atoms should 
be maximal and one needs to avoid hypercoordination of the C or O atoms. In CLi4 or OLi2 all 
valence electrons are involved in C— Li and O— Li bonds, respectively, maximizing the charge on 
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Figure 1: Optimized geometries of (a) CLi4H8, CLi4H24, and (b) OLi2H4, OLi4H2o complexes. 
Carbon, oxygen, lithium and hydrogen atoms are shown as black, red, yellow and blue spheres, 
respectively. 

the Li atoms. The most stable geometries of CLi4 and OLi2 have a high symmetry, i.e. and Dooh, 
respectively. All Li atoms in a molecule are then equivalent. 

A Bader charge analysis of the calculated self-consistent charge densities gives a charge of 
+0.82^ on each of the Li atoms in CLi4, and a charge of +0.88^ on the Li atoms in OLi2.^^ These 
high charges facilitate the electrostatic bonding of hydrogen molecules to the Li atoms in these 
molecules, similar to that in isolated Li+ ions.^^ Hydrogen molecules then cluster around the Li 
sites of the CLi4 and OLi2 molecules. Examples of the structures of hydrogenated molecules are 



shown in Figure 1 and bond lengths before and after hydrogenation are given in Table 1 



Table 1: Calculated bond lengths (A), Bader charges (e), hydrogen binding energies E^^ (eV/H2) 
and gravimetric hydrogen densities of CLi4 and OLi2 molecules. For Li— H2 the bonding distance 
is measured between Li atom and the center of H— H bonds. Average numbers are given and the 

numbe rs between () indicate the spreads. 

bond lengths (A) 



Molecule 


C,0-Li 


Li-H2 


H-H 


charge on Li 




H wt. % 


CU4 


1.898(3) 






+0.819(2) 






CLi4H8 


1.888(5) 


2.036(12) 


0.756(0) 


+0.827(2) 


0.14 


16.86 


CLi4H24 


1.908(3) 


2.158(7) 


0.763(2) 


+0.851(3) 


0.10 


37.82 




1.643(0) 






+0.876(0) 






OLi2H4 


1.671(1) 


1.862(5) 


0.793(0) 


+0.869(0) 


0.21 


11.89 


OLi2H20 


1.711(1) 


2.140(144) 


0.764(13) 


+0.881(1) 


0.13 


40.29 
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Also given in [Table 1| are the hydrogen binding energies, which are calculated by subtracting 
the total energy of the hydrogenated polylithiated molecules from the sum of the total energies 
of the isolated polylithiated molecules and the hydrogen molecules, divided by the number of 
hydrogen molecules. Each Li atom in a CLi4 molecule can bind three H2 molecules as shown in 



Figure 1[ which leads to 37.8 wt.% H. The average binding energy decreases from 0.14 eV for the 
first H2 molecule (per Li atom), to 0.10 eV for the fully hydrogenated state. The decrease reflects 
a net repulsive interaction between the H2 molecules surrounding a Li atom. Steric hindrance 
prevents binding a fourth H2 (at least with a binding energy > 0.04 eV). The H2 binding energies 
are consistent with those found for isolated Li^ ions, bearing in mind that the Li atoms in CLi4 
have a somewhat reduced charge.^ 

OLi2 molecules bind hydrogen more strongly than CLi4 molecules. As the net charge on a Li 
atom in OLi2 is larger than in CLi4 this can be expected if the binding is electrostatic. Because 
the Li atoms in OLi2 are less crowded than in CLi4, each Li atom can accommodate up to five H2 
molecules, before steric hindrance prevents binding of additional hydrogen. The average binding 
energy decreases from 0.21 eV for the first H2 molecule (per Li atom) to 0.13 eV in the fully 
hydrogenated state. It leads to 40.3 wt.% H for the OLi2H2o complex. 

The fully hydrogenated CLi4H24 and OLi2H2o complexes keep a high symmetry, i.e. and 



D2d, respectively. In intermediate hydrogenation stages the symmetry can be broken, see [Figured 
After hydrogenation the C— Li and O— Li bonds are stretched by 1% respectively 4%. This is ac- 
companied by a slight increase of the partial charge on Li by 4% respectively 0.6%. The optimized 



calculated bond length of an isolated H2 molecule is 0.748 A. [Table 1| shows that in the complexes 
the H— H bond lengths increase by 1-6%. The limited size of these changes agrees with the no- 
tion of a weak interaction with the H2 molecules. The absolute size of these changes generally 
correlates with the binding energy. 
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Figure 2: Optimized geometries of (a) (CLi3)2, (b) (CLi4)2, (c) (OLi)2 and (d) (OLi2)2. 

Adsorption on Graphene 

It is clear that individual CLi4 and OLi2 molecules can bind a substantial number of H2 molecules. 
In order to maintain reversibility during repeated hydrogen loading and unloading cycles, one 
should however prevent the formation of clusters or aggregates of CLi4 or OLi2 molecules. From 
a thermodynamic point of view such clusters could form. For instance, the equilibrium geom- 



etry of C2Li8 is shown in Figure 2 ^b).^^ The energy of this structure is 2.54 eV/CLi4 lower in 
energy than that of two separate CLi4 molecules. Similarly, two OLi2 molecules can form a clus- 
ter with a binding energy of 1.37 eV/OLi2, see [Figure^ d). Presumably also structures of larger 
clusters exist that have a favorable binding energy. One could rely on kinetic barriers preventing 
clustering of polylithiated molecules. Here we consider alternative possibilities, which consist of 
immobilizing the molecules by binding them to a holder template. Microporous carbon and carbon 
nanostructures would be obvious substrates for these molecules. 

As a prototype system we consider binding to a single sheet of graphite, i.e. to graphene. We 
use a 4x4 graphene supercell, which contains 32 atoms and has an optimized cell parameter of 
9.866 A. This cell is sufficiently large to prevent a direct interaction between molecules bonded to 
graphene in neighboring cells. The calculated binding energies of single CLi4 and OLi2 molecules 
bonded to graphene are 0.73 and 0.55 eV, respectively. Both of these energies are considerably 
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Figure 3: Top and side views of the optimized geometries of graphene binding (a) two CLis 
molecules, or (b) four OLi molecules per supercell, in their fully hydrogenated forms. 

smaller than the corresponding molecular dimerization energies, indicating that the interaction 
with the substrate is not sufficiently strong. One has to consider increasing this interaction. 

For this we pursue two different strategies consisting of modifying the molecules or the sub- 
strate, respectively. Removing one Li atom from CLi4 or OLi2 strengthens the bonding to graphene. 
The calculated binding energy between CLis and graphene is 2.00 eV. CLis molecules on graphene 
are preferably positioned over the center of C— C bonds. Both faces of graphene can be function- 
alized, which leads to a geometry as shown in [Figure ^ a). This structure is stable, but it is not 



the absolute minimum energy structure. C2Li6 in its equilibrium structure, cf. Figure 2 ^a),^^ has 
a binding energy of 3.86 eV/CLis with respect to two separate CLis molecules. It means that, 
although the CLis— graphene structure is kinetically stable, there is still a thermodynamic driving 
force towards aggregation of CLis molecules. 

OLi molecules also bind strongly to a graphene substrate. Since an OLi molecule is smaller 
than a CLis molecule, it is possible to quadruple the packing density of these molecules, as is 



shown in [Figure 3[ b). OLi molecules preferentially bind to the substrate through the formation 
of an O— C bond. The calculated binding energy is 2.22 eV, which is comparable to the 2.26 
eV/OLi aggregation energy of two OLi molecules into 02Li2. The 02Li2 molecule has a planar 



structure as shown in Figure 2 c). This means that the OLi-graphene structure is close to being 
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Figure 4: Top and side views of the optimized geometries of the (a) Be doped and (b) B doped 
graphene structures. 

thermodynamically stable. 



Adsorption on Doped Graphene 

Another possible way of increasing the binding between polylithiated molecules and the substrate 
is by chemical modification of the graphene layer. Introducing an impurity in graphene might 
for instance not only strengthen the interaction with CLi^ and OLi^ molecules, but could also 
maintain the positive charge on the Li atoms, and therefore the binding to H2 molecules. Here we 
consider modification of graphene by atomic substitution of C by B or Be atoms. Previous studies 
have shown that B doping of graphene does not cause any major deformation of the graphene 
structure.!^ B doping concentrations of as high as 17 at. % can be achieved in thin films.l^From 
experiment it has been suggested that the preferential geometry upon high B doping has two B 



atoms in para positions of an hexagonal ring in the graphene plane, as shown in [Figure 4 ^b) 



Indeed, the stability of this geometry has recently been confirmed by a theoretical study Here 
we consider such a double B substitution in our 4x4 graphene supercell, which gives a doping 
concentration of 6.25 at. %. 

As an alternative to B doping we also consider substituting C with Be atoms. Replacing one 
C with a Be atom does not alter the planar structure of graphene. It is energetically advantageous 
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Figure 5: Side views of the fully hydrogenated structures of CLi4 on (a) Be doped and (b) B doped 
graphene and of CLi5 on (c) Be doped and (d) B doped graphene. 

however to replace a single C atom by a pair of Be atoms. The two Be atoms are then displaced 



above and below the graphene plane by 0.855 A, as shown in [Figure 4[ a), whereas the rest of 
graphene maintains its original structure. The Be atoms can act as pedestals that increase the 
distance between the graphene plane and the molecules we are going to attach, which could limit 
steric hindrance. All of the Be or B doped graphene structures have a nonmagnetic ground state. 

Modifying graphene by Be or B doping increases its interaction with polylithiated molecules. 
The molecules preferentially bind to the dopant atoms and form C,0— B,Be bonds. From the point 
of view of binding strength, there is no need to modify the molecules by extracting a Li atom 
before attaching them to the substrate. CLi4 binds to a B or a Be site with binding energies of 2.91, 
3.81 eV, respectively. The binding geometries are shown in [Figure 5| The binding between OLi2 
and B or Be doped graphene has an energy of 2.49, 3.76 eV , respectively. Since all these binding 
energies are well above the molecular dimerization energies one expects that aggregation of these 
molecules is now suppressed. Remarkably, doped graphene also forms strong bonds to other CLi^ 
and OLim molecules, including hypercoordinated ones such as CLi5, OLis and OLi4, leading to 



similar binding geometries. The binding energies are given in [Table 2| and examples of binding 
geometries are shown in [Figure 5 



With the polylithiated molecules strongly bonded to (doped) graphene we now study the hy- 
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CLi3 


CLi4 


CLis 


OLi 


0Li2 


OLi3 


OU4 


undoped 


2.00 


0.73 




2.22 


0.55 






Be doped 


4.12 


3.81 


3.96 


4.14 


3.76 


3.55 


3.00 


B doped 


4.05 


2.91 


3.23 


3.45 


2.49 


3.13 


3.69 



drogenation properties of these systems. For CLis and OLi bonded to undoped graphene we find 
that each of the Li atoms in the molecules above and below the graphene plane can bind three H2 
molecules, see [Table 3| The average hydrogen binding energy is 0.10 and 0.12 eV, respectively, 



which is comparable to that in free CLi4 and OLi2 molecules, compare Table 1[ The resulting ge 



ometries are given in [Figure 3| The gravimetric densities in the CLis- and OLi-graphene systems 
are 7.46 and 7.85 wt.% H, respectively. 

As discussed above, doping of graphene with Be or B increases the binding of polylithiated 
molecules to the substrate. The effect of doping on the binding energy between hydrogen and 
these molecules is small in most cases, as can be observed in [Table 3[ Only for H2 on CLis and 



OLi bonded to Be doped graphene the binding energy is significantly increased. The distance 
between the hydrogen molecules and the graphene plane is larger in these systems, as compared 
to the other cases. The gravimetric hydrogen density in the doped graphene systems is determined 
by the concentration of dopant atoms, since the polylithiated molecules are bonded to these atoms, 
as well as by the number of hydrogen molecules bonded to a single polylithiated molecule. 

The latter is foremost determined by steric hindrance. A Li atom that is close to the graphene 
plane or to other Li atoms can in general capture only one or two H2 molecules, whereas the ones 
that are further away from the plane and other Li atoms can bind three H2 molecules. This is 



illustrated in [Figure 5[ We have fixed the doping concentration of graphene. In particular for the 
OLi^ systems one can still increase the doping concentration without introducing too much steric 
hindrance, which would increase the gravimetric density substantially. 

Hydrogen binding energies decrease slightly for hypercoordinated polylithiated molecules, i.e. 
CLi^, n > 3 and OLi^, m > 1. [Figure 5[ shows that this leads to a trade-off between hydrogen 
binding energy and gravimetric density. Nevertheless it is clear that there are several possible 
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Table 3: Average hydrogen binding energies Ej^ ^ (eV/H2), number of hydrogen molecules #H2 
bonded per CLi^ or OLi^ species, and gravimetric hydrogen densities of polylithiated molecules 
bonded t o (doped) graphene. 



undoped Be doped B doped 







#H2 


H wt.% 




#H2 


H wt.% 


'^b 


#H2 


H wt.% 


CLi3 


0.10 


9 


7.46 


0.15 


7 


5.83 


0.11 


8 


6.72 


CLi4 








0.09 


9 


7.17 


0.09 


9 


7.29 


CLis 








0.08 


11 


8.40 


0.07 


11 


8.53 


OLi 


0.12 


3 


7.85 


0.23 


3 


2.70 


0.10 


3 


2.75 


0Li2 








0.10 


6 


5.10 


0.09 


6 


5.19 


0Li3 








0.12 


6 


4.96 


0.10 


6 


5.04 


0Li4 








0.08 


9 


7.06 


0.06 


9 


7.18 



systems and structures that attain a gravimetric density higher than 5% H with H2 binding energies 
in excess of 0.1 eV. 



Conclusions 

In conclusion, we have shown by means of first-principles DFT calculations that polylithiated 
molecules can be used as versatile building blocks for hydrogen storage materials. Li— C and 
Li— O bonds are strongly polar and the partially charged Li atoms can attract up to three to five 
H2 molecules with an average binding energy in the range 0.10-0.15 eV. CLi4 and 0Li2 molecules 
can then bind 12 and 10 hydrogen molecules, respectively, giving very high gravimetric densities 
of 37.8 and 40.3 weight % H. Binding polylithiated molecules to carbon nanostructures such as 
graphene can be used to immobilize them. In particular doping graphene with B or Be increases 
the binding energy. Although there is a trade-off between the size of the H2 adsorption energy 
and its gravimetric density in these systems, densities in the range 5-8.5 wt.% H can be reached 
without seriously decreasing the adsorption energy. 
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